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ABSTRACT 

We measure the correlation between sky coordinates of the Swift BAT catalogue of active 
galactic nuclei with the arrival directions of the highest energy cosmic rays detected by the 
Auger Observatory. The statistically complete, hard X-ray catalogue helps to distinguish be- 
tween AGN and other source candidates that follow the distribution of local large-scale struc- 
ture. The positions of the full catalogue are marginally uncorrelated with the cosmic ray arrival 
directions, but when weighted by their hard X-ray flux, AGN within 100 Mpc are correlated 
at a significance level of 98 per cent. This correlation sharply decreases for sources beyond 
^ 100 Mpc, suggestive of a GZK suppression. We discuss the implications for determining 
the mechanism that accelerates particles to these extreme energies in excess of 10^^ eV. 
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1 INTRODUCTION 

The recent announcement of anisotropy in the arrival directions of 
ultra-high energy cosmic rays (UHECRs) by the Auger Collabora- 
tion (Abraham et al. 2007, 2008) has yielded much interest. By cor- 
relating the UHECR arrival directions with the positons of nearby 
active galactic nuclei (AGN) on the sky, they suggest that AGN are 
responsible for accelerating these particles. 

Modern experiments are approaching the size and sophis- 
tication needed to measure the energy and direction of UHE- 
CRs with enough precision to reliably determine their sources. 
The Auger Observatory combines measurements of Cerenkov ra- 
diation from particle interactions in surface detectors and flu- 
orescence from molecu l es in the atmosphere excited by the 
cascade j Abraham et alj 2004 ). For r e views o f earlie r exper - 
imental work, see iNagano & WatsonI ilOOd) : ICronid ( l2005h : 
ISokolskv & ThomsonI ^200% ! 

The origin of particles with energies > 10^^ eV has been a 
longstanding mystery in high energy astrophysics (see Hillas 1984 
for a review). Supemovae, gamma-ray bursts, pulsars, shock fronts 
in galaxy clusters, decays of exotic massive particles, and various 
classes of AGN have all been proposed as acceleration or gen- 
eration sites. All but the highest energy cosmic rays are thought 
to be significantly deflected by Galactic magnetic fields, which 
makes it difficult to trace arrival directions back to astronomical 
sources. Numerous attempts have been made to correlate differ- 
ent source catalogues with events observed by cosmic ray detec- 
tors including Yakutsk, Fly's Eye, AGASA, HiRes, and now Auger 
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(e.g.. IStanev et alj Il995l: Isingh et alj |2004 [Hague et alJ l2007l : 
iGorbunov et alj2007l :'Kashti & Waxmanl20 08":' Ivano v et al.l2008h . 
Related studies have analysed the auto-correlation of sources and 
events as well as the occurrence of multiple co-located UHECR 
detections to place constraints on the number of sources responsi- 
ble for these events (Takeda et al. 1999; C uoco et al. 2008) . AGN 
have long been known to be capable of gen erating electric pot ential 
differences in excess of 10^^ volts ( Bums & Lovelacj|l982h , and 
given the further suggestion from the Auger team that their posi- 
tions are correlated with UHECR arrival directions, we focus only 
on these objects. 

The correlation of Auger events with extragalactic sources 
given bv'Abraham et al.l j2008l) is based on an inhomogeneous list 
of AGN collected from the literature ( Veron-Cettv & Veron 20061, 
hereafter V-C). As pointed out by others dKashti & Waxmanl l2008h 
and explained in the companion paper to the catalogue itself, the 
compilation is not complete and should not be used for statistical 
purposes. While using such a list can be adequate for disproving the 
hypothesis of an isotropic distribution of arrival directions, positive 
proof of correlation with sources calls for a statistically complete 
sample. Additionally, the distribution of local AGN is similar to that 
for all local large-scale structure. A complete sample of sources se- 
lected for a certain property can distinguish between different trac- 
ers of large-scale structure and also determine if a particular factor 
is relevant to UHECR acceleration. 

In this letter, we study the correlation of Auger UHECR arrival 
directions with the po sitions of AGN fro m the hard X-ray selected 
Swift BAT catalogue ffueller et al.ll2008l) . Selection from the hard 
X-ray band reduces the bias due to absorption that impacts lower 
energy bands. Optically selected lists, including many that are com- 
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piled into the V-C catalogue, are less likely to include the obscured 
nuclei that make up an increasing fraction of the AGN popula tion 
at low redshifts ( iBrandt & Hasingedl2005l) . iMatt et alJ kood) de- 
duce from hard X-ray observations that heavily obscured sources 
must locally outnumber unobscured sources by a large factor; the 
nearest 3 AGN are all highly obscured. 

Since hard X-ray observations are less sensitive to obscura- 
tion, they provide a more accurate indicator of the accretion rate 
and the intrinsic luminosity of an AGN. We aim to test the idea that 
the intrinsic luminosity of an AGN is a sign of its particle acceler- 
ating ability by measuring the relation between hard X-ray flux and 
the positional correlation of AGN with UHECR arrival directions. 
In the next section, we will present the catalogues of UHECRs and 
AGN. Sj3]will describe the statistical tests used to measure corre- 
lation between the two lists as well as the results of these tests. 
We discuss the implications of these results and compare them to 
previous analyses in Sj?) and ij5]concludes the paper. 



2 DATA 

In this section we discuss the catalogues of UHECRs and AGN 
used for testing the correlation between events and sources. We do 
not combine data from different cosmic ray observatories because 
of possible inconsistencies in energy calibration and angular reso- 
lution. A statistically complete, hard X-ray selected sample of AGN 
will reduce the bias toward unabsorbed sources present in the V-C 
catalogue, which has mainly been selected from optical observa- 
tions. This sample should clarify whether the accretion processes in 
AGN that generate hard X-rays are relevant to high energy particle 
acceleration. By narrowing the sample from a large and inhomoge- 
neous list of AGN, we can test the possibility that the observed cor- 
relation of UHECRs and AGN is merely due to the fact that AGN 
are distributed similarly to other sources in the local structure. 

2.1 UHECR events 

The Auger Observatory ^Abraham et alj|2004l2008h has been op- 
erating stably in Argentina since 2004, using a hybrid system of 
telescopes to measure fluorescence in the atmosphere and water 
tanks to detect Cerenkov light from relativistic particles. The avail- 
able UHECR list is comprised of 27 events with energies above 
5.7 X 10^^ eV from an integrated exposure of 9000 km^ sryr. 

The relative exposure is independent of energy in this range, 
nearly u niform in right asc ension, and has a declination dependence 
given bv lSommerjhoOlh . The latitude of the Auger Observatory 
is —35.2° and it has a maximum zenith angle acceptance of 60°. 

The event arrival directions are determined with an angular 
resolution of better than 1°. However, magnetic fields of unknown 
strength will deflect charged particles on their trajectories through 
space. The advantage of studying the highest energy events is that 
this deflection is mi nimized, but it can still be up to ~ 10° in the 
Galactic field (e.g., IStanevlll997h . Magnetohydrodynamical simu- 
lations of extragalactic fields have produced conflicting predictions 
for U HECR deflections, rangin g from tens of d egrees cSigl et al.l 
|2004) to less than a few degrees jPolag et all2005l) . 

2.2 AGN catalogue 

We compare the arrival directions of Auger UHECRs with the lo- 
cations of AGN in the SWIFT BAT survey. Details of the survey 



and result s from the first 3 and 9 mo nths of observations are pre- 
sented by Markwardt et alj ( l2005h and lTueller et al.l ( l2008h : we use 
a catalogue updated with the first 22 months of data. 

The list includes 254 objects that have known redshifts and 
signal to noise ratios greater than 4.8, corresponding to one false 
detection across the sky. Source identification is incomplete at 
low Galactic latitudes (|&| < 15°), but the hard X-ray cover- 
age (14 — 195 keV) provides a uniquely unbiased sample of lo- 
cal AGN over the rest of the sky. The median astrometric un- 
certainty is 1.7 arcmin and the limiting sensitivity is a few times 
10^^^ ergcm"^ s~^. All but four sources in the 22 month cata- 
logue have X-ray positions and spectra from previous observations 
or follow-up with the SWIFT XRT 

We restrict our statistical analysis to the complete sample 
within the Auger field of view (S/N > 4.8, |6| > 15°, (5 < 24.8°), 
leaving 138 AGN in the catalogue. For comparison with the orig- 
inal correlation between Auger events and AGN, there are 57 
Swift sources within the Auger field of view at distances less than 
100 Mpc, and 6 of these are not in the V-C catalogue. The V-C cat- 
alogue has 410 objects within 100 Mpc and the Auger field of view 
outside of the Galactic plane. Throughout this paper, distances are 
determined using a cosmology of Ho = 71 kms^^ Mpc^^ and 
fim = 0.27 with a flat universe. 



3 ANALYSIS 

Figure[T]shows the locations of Auger UHECRs and AGN from the 
BAT catalogue in supergalactic coordinates. The relative exposure 
of Auger observations is plotted in yellow contours, and the AGN 
are shaded blue according to the product of their X-ray flux and 
the relative Auger observation time. Blue curves indicate where the 
BAT catalogue is incomplete at low Galactic latitudes, j6j < 15°. 

Visually, we can see several interesting correlations and non- 
correlations. Cen A is the brightest source and has two Auger 
events within 3°. The likelihood of observing 2 of the 27 events 
within this angular separation of the brightest AGN in the entire cat- 
alogue for an isotropic distribution of arrival directions in Auger's 
field of view is about 3 x 10"*. The 6 AGN within 20 Mpc are 
marked with white crosses, and 4 of these (including Cen A) have 
UHECR arrival directions within 6° . If this proportion of observed 
UHECRs does in fact arrive from nearby AGN and the sources are 
roughly isotropic on large scales, we can infer that there must be 
significant propagation losses for cosmic rays from distant sources. 
However, clusters of 2 and 3 Auger events have no corresponding 
AGN, and the overall correlation is unclear. Certainly there is not a 
1-to-l correspondence between all of the observed UHECRs and a 
hard X-ray emitting AGN. 

To quantify the likelihood that the two sets of locations 
arise from the same population of source s, we employ the two- 
dimensional generalizat ion of the K-S test iFasano & Franceschinil 
ll987l ; |Press^ai]|2007h . This test can be used to compare a data 
sample against a model distribution, or to compare the properties 
of two data samples against each other. To test the null hypothe- 
sis that a given class of sources is responsible for UHECRs, one 
could construct a model of the expected cosmic ray flux map using 
a catalogue of positions, intrinsic luminosities, distances, expected 
absorption losses, and the detector acceptance as a function of po- 
sition. Additionally, one should account for interactions with back- 
ground radi ation which produce propagation losses known as the 
GZK effect l lGreiserJl96^ : IZatsepin & Kuz'minll9"66h . lSingh et al.l 
^2004.) used this approach to measure the correlation of infrared 
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Figur e 1. Map of Auger UHECRs (open red circles) and BAT AGN within 100 Mpc (filled blue circles) in supergalactic coordinates jde Vaucouleurs et alj 
Il976h . The blue colour depth is scaled by the hard X-ray flux and Auger exposure, relative to Cen A. The 6 AGN in the catalogue within 20 Mpc are marked 
with white crosses, with Cen A at (159.7°, —5.2°). Yellow contours have equal integrated exposures. Blue boundaries show where the AGN catalgoue is 
incompete due to the Galactic plane, \b\ < 15°. 



galaxies with UHECR events observed by AGAS A, including mod- 
els for the injection spectra and propagation effects. 

We might expect the hard X-ray flux, as a measure of the in- 
trinsic power of the source, to be a reasonable proxy for the cosmic 
ray flux. But because a model of GZK effects requires an input 
spectrum of cosmic rays and the injection mechanism is unknown, 
we opt to begin by simply comparing the AGN locations with the 
UHECR arrival directions. Subsequently, we weight these positions 
by the product of the X-ray flux and the relative Auger exposure. 

The two-dimensional K-S test measures the fraction of data 
points lying in the natural quadrants defined around each data 
point. The statistic D is the maximum difference between the frac- 
tions of each of the two data sets that lie in a quadrant, found 
after iterating over all points and their respective quadrants. The 
strength of the correlation between two catalogues is reported as 
the integral probability distribution P{D^/n > observed), where 
n = NiN2/{Ni + A''2), and TVi and N2 are the number of data 
points in the two sets. This measurement can be used to deter- 
mine the similarity of sets of positions on the sky. Statisticafly 
speaking, high values of P do not prove the null hypothesis that 
the two catalogues arise from the same population, but low values 
{P < 20 per cent) call for the rejection of that null hypothesis. 

The probability that the two sets of data are from the same 
population can be determined with an analytical approximation if 
the number of data points is sufficient, but due to the limited num- 
ber of events available and the irregular shape of the allowed coor- 
dinate plane, we generate Monte Carlo simulations of cosmic ray 
arrival directions to compare with the observed events. For each set 
of AGN considered, 10,000 lists of positions are randomly gener- 
ated in a distribution on the sky that is proportional to the Auger 
exposure, excluding low galactic latitudes. Each list has the same 
number of UHECRs as observed in the selected region, and has its 



D statistic calculated in the same way as the real data. The sig- 
nificance of the correlation between a catalogue of AGN and the 
observed UHECRs is reported as P[D^Jn > observed), which 
is the percentage of randomly generated UHECR lists that have a 
higher value of D^/n than the real data under consideration. Thus, 
extreme values of P suggest that the observed events are not dis- 
tributed isotropically, and high (low) values indicate a good (poor) 
correlation between the Auger UHECRs and the given AGN list. 

It is important to ensure that the test does reliably differentiate 
between a set of points correlated with a reference catalogue and 
an isotropically distributed set. We picked random subsets of AGN 
from the larger V-C catalogue to serve both as a list of AGN and 
a list of UHECRs, with the number of sources and events equal to 
those in the actual tests on Swift and Auger data. Average values for 
the correlation showed f > 70 per cent, setting a rough threshold 
for accepting that two data sets belong to the same distribution. 

The test differs from the one performed by the Auger team in 
that the angular distance between the source position and arrival 
direction is not a parameter of the fit, so it is not as sensitive to the 
undetermined size of magnetic deflections. The angular scale of the 
test's sensitivity is not clearly defined and depends on the separa- 
tion of points in the data set, but Monte Carlo simulations indicate 
that correlations decrease monotonically with angular separation, 
and that this decline has a scale length of approximately 10 — 15°. 
Additionally, we use the full energy range of the published Auger 
UHECR catalogue, which has a threshold selected to maximize the 
correlation in their analysis. We allow source distances in the full 
range of the AGN catalogue, though completeness becomes an is- 
sue at higher distances. 

In our application, we use the equatorial positions for the 
Auger events and Swift AGN as the data points. Restricting the data 
sets to regions and significance levels where the BAT catalogue is 
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Figure 2. Redshift dependence of the correlation between AGN positons 
and UHECR arrival directions. Light (dark) boxes show the unweighted 
(flux-weighted) values of P, defined in the text. In the most distant bin, 
flux-weighting decreases the significance of the correlation. There are 30 
AGN in each comoving distance bin. 



complete > 15°,S/N 4.8) and to where the Auger exposure 
is non-zero {5 < 24.8°) leaves 138 AGN and 19 UHECRs. For 
these catalogues, we measure P = 50 per cent. This value jumps 
to 98 per cent after weighting the AGN coordinates by their hard X- 
ray flux and relative Auger exposure. If we cut the AGN catalogue 
to those 57 with distances less than 100 Mpc, the unweighted prob- 
ability is 84 per cent. The flux-weighted value is unchanged, indi- 
cating that the correlation is dominated by nearby, bright sources. 
However, the removal of Cen A does not alter the correlation val- 
ues by more than a few per cent. Cuts on the catalogue to dis- 
tances smaller than 100 Mpc do increase the correlation slightly, 
but quickly the number of sources becomes too small to obtain re- 
liable values from the test. For comparison, the 410 sources in the 
V-C catalogue with d < 100 Mpc, > 15°, and S < 24.8° 
produce a correlation with the UHECRs of only 55 per cent. 

To determine the redshift dependence of the correlation, we 
test 4 redshift-sorted bins of 30 AGN each, with mean comoving 
distances of 35, 80, 131, and 248 Mpc. The results are plotted in 
Figure |2l which shows a clear decrease in the correlation at larger 
distances. Replicating the plot with different bin sizes, we consis- 
tently see a sharp decrease in the correlation near 100 Mpc. 

We performed a similar test on the luminosity-dependence 
of the P values and found that the correlation in the full cata- 
logue is actually lower for more luminous sources than for dimmer 
ones. It is possible that the strong radiation field surrounding lu- 
minous AGN interferes with UHECR propagation (Dermer 2007 ). 
But when restricting the AGN to distances less than 100 Mpc, the 
trend reverses and more luminous sources tend to correlate better 
than less luminous ones. In both cases, the most powerful AGN 
are at a higher average distance than dim ones. We interpret these 
results to imply that the cosmic ray flux is correlated with the lu- 
minosity of AGN, but propagation losses on a distance scale of 
~ 100 Mpc interfere with the effect. 

We would also like to test the correlation for subclasses of 
AGN. There are 46 "Type 1" objects (Seyfert 1/1.2 galaxies) and 
53 "Type 2" objects (Seyfert 1.8/1.9/2 galaxies) at high Galactic 
latitudes and within the Auger exposure. There are not enough ob- 
jects in other subclasses such as blazars or quasars to obtain statis- 
tically significant results. For Type 1 objects, the unweighted and 
weighted P values are 59 and 96 per cent, respectively, while for 
Type 2 objects these values are 53 and 99 per cent. For both types, 
excluding sources at distances greater than 100 Mpc increases the 



unweighted correlation by --^30 percentage points, while the flux- 
weighted values remain the same. 

The overall results are consistent with visual expectations, 
in the sense that the correlation between the full catalogues is 
marginal at best, but various cuts and weighting factors can im- 
prove the correspondence. We risk falling into the trap of using 
a posteriori cuts to find the best matches and reporting prejudical 
results. Lengthy arguments have taken place over the usefulness of 
applying "penalty fa ctors" to account for th e statistical biases pro- 
duced by these cuts (Evans et al .120031 ; iFirUev & Westerhofl 2004]) . 
but we omit this type of analysis and merely point out that the stated 
probabilities should be taken with these caveats in mind. 



4 DISCUSSION 

Our results from the K-S tests do not cause us to reject the null 
hypothesis that AGN are the source of UHE CRs, in disagreement 
with the conclusion of lGorbunov et al.l ( |2007|) . On the contrary, we 
find significant correlations that sugge st that AGN are the so urces 
of UHECRs, supporting the claims of [Abraham et alj ( l2007h . The 
hard X-ray catalogue used in this paper produces a higher correla- 
tion with UHECR arrival directions than the larger V-C catalogue 
which was previously tested. From this information alone, it is still 
possible that another type of source with a local distribution similar 
to that of AGN is responsible for UHECR acceleration. But when 
weighting by the hard X-ray flux or selecting the most luminous 
nearby sources, the measured correlation between AGN and UHE- 
CRs increases. Thus, a confounding class of sources would have 
to trace both the local spatial and luminosity distributions of AGN, 
which seems unlikely. 

If AGN are positively identified as the sources of UHECRs, 
the second main result of this paper is the significant decline in cor- 
relation at a distance of ~ 100 Mpc. This drop is predicted from 
the GZK cutoff, which suggests that the flux of cosmic rays in the 
energy range considered here is depleted due to photopion produc- 
tion and pair creation from interactions with the cosmic microwave 
background. For cosmic rays with energies above 5 x 10^^ eV, the 
observable "horizon" is expected to be of order 100 Mpc, with de- 
viations due t o the dis tribution of sources and the masses of accel- 
erated nuclei jHarari e t al. 2006). HiRes has detected a suppression 
of the cosmic ray spectrum at ener gies above 6 x 10^^ eV which has 
been attributed to the GZK cutoff jAbbasi et alj2008l) . Though this 
detection has been supported by the Auger data, it is conceivable 
that the intrinsic spectrum of cosmic rays falls off at this energy due 
to some cutoff in the acceleration mechanism. An observation that 
UHECRs arrive from nearby sources and not distant ones would 
provide direct evidence for the GZK effect. 

The drop observed in the correlation between AGN positions 
and UHECR arrival directions at 100 Mpc, as shown in Figure|2l is 
suggestive of such a suppression. It may arise due to other effects 
including increased angular deflections due to extragalactic mag- 
netic fields and incompleteness of the AGN catalogue beyond this 
distance (Abraham et al. 2008). However, if magnetic deflections 
are the culprit we would expect the correlation to begin decreasing 
before 100 Mpc, which is not observed. Catalogue incompleteness 
affects the detection of less luminous AGN at large distances, but 
for local AGN the correlation reduces with faintness. Additionally, 
the evidence discussed earlier that 4 nearby AGN have associated 
events does, from Olber's paradox considerations, strongly argue 
that the sources of all detected UHECRs are at very low redshift. 

The association of two UHECRs with the brightest AGN 
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in the Swift catalogue allows for speculation about the nature 
of particle acceleration mechanisms. Several studies have sug- 
gested that blazars are lik ely candidates for UHECR sources (e.g., 
iTinvakov & Tk achev 2001), though the Auger data do not support 
a correlation (jHarari et al. 2007). Cen A is a nearby jetted AGN 
with two large radio lobes visible from the side. If Cen A is in fact 
responsible for some of the observed UHECRs, then cosmic rays 
do not need to be emitted directly along a jet toward the observer. It 
is conceivable that jets may still be important for cosmic ray accel- 
eration and that the events seen towards Cen A were deflected by 
strong magnetic fields surrounding the AGN. Our division of the 
AGN catalogue into "Type 1" and "Type 2" classes yielded similar 
results from the K-S test for correlation, but further tests on dis- 
tinct classes of AGN, including more heavily obscured sourcefl 
will help clarify which mechanism or orientation is necessary to 
produce the observed UHECRs. 

As data for a larger number of UHECRs become available, we 
will be able to further differentiate between the properties of ac- 
celeration sites. We have shown that the correlation between AGN 
positions and UHECR arrival directions does vary with the hard 
X-ray flux and the distance to the source. It will be interesting to 
test other properties to constrain models for the acceleration mech- 
anism. Black hole spin may be an important parameter but has not 
yet been tested because of the difficulty in measuring its value for 
a sizable sample of sources. Man y models exist for the extra ction 
of energy from black holes (e.g., iBlandford & Znai^ll977l) . and 
a further analysis of the influence of the model parameters on the 
correlation between source positions and UHECR directions will 
more precisely determine which sources are accelerating these par- 
ticles. Conversely, when the cosmic ray sources are clearly identi- 
fied, more detailed observations in the realm of particle astronomy 
will improve our understanding of the physical mechanisms taking 
place in these sources. 



5 CONCLUSION 

We have measured the correlation between the positions of the 
Swift hard X-ray selected catalogue of AGN and the arrival direc- 
tions of UHECRs observed by Auger. No significant correlation is 
found for the full catalogue, but imposing weights based on the ex- 
pected cosmic ray flux results in a correlation with a significance 
level of 98 per cent. A steep drop in the correlation for sources 
beyond 100 Mpc suggests that the GZK effect does suppress the 
cosmic ray flux at these energy and distance scales. 

Though our study lends support to the claims of the Auger 
Collaboration that AGN are responsible for UHECRs, the connec- 
tion should still be tested with better statistics as observations con- 
tinue. The identification of 2 UHECRs from the direction of Cen A 
and the improvement of the correlation when weighting by hard X- 
ray flux offer tantalizing hints about the acceleration mechanism. 
Though the nature of these exotic UHECRs remains largely un- 
known, future studies in particle astronomy may lend insight into 
other areas such as Galactic and extragalactic magnetic fields and 
the physics behind AGN engines. 



We searched for a correlation with several heavily obsc ured sources too 
faint t o be detected in the BAT catalogue, selected from iLevenson et al] 
We have not identified any obscured AGN located near the UHECR 
arrival directions, but the sample size is too small to draw conclusions about 
this population of AGN in relation to the UHECRs. 
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